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ABSTRACT 



We develop a method for reconstructing the two-dimensional, projected gravitational potential of galaxy clusters from observed line- 
of-sight velocity dispersions of cluster galaxies. It is the third of an intended series of papers aiming at a unique reconstruction method 
for cluster potentials combining lensing. X-ray, Sunyaev-Zel'dovich and kinematic data. The observed galaxy velocity dispersions 
are deprojected using the Richardson-Lucy algorithm. The obtained radial velocity dispersions are then related to the gravitational 
potential by using the tested assumption of a polytropic relation between the effective galaxy pressure and the density. Once the 
gravitational potential is obtained in three dimensions, projection along the line-of-sight yields the two-dimensional potential. For 
simplicity we adopt spherical symmetry and a known profile for the anisotropy parameter of the galaxy velocity dispersions. We test 
the method with a numerically simulated galaxy cluster and galaxies identified therein. We extract a projected velocity-dispersion 
profile from the simulated cluster and pass it through our algorithm, showing that the deviation between the true and the reconstructed 
gravitational potential is < 10 % within « l.lh^'Mpc from the cluster centre. 

Key words. (Cosmology:) dark matter. Galaxies: clusters: general. Galaxies: clusters: kinematics. Gravitational lensing: strong. 
Gravitational lensing: weak 



1. Introduction 

The success of the cosmological standard model implies that 
massive, gravitationally bound cosmic objects such as galaxy 
clusters should be dominated by dark matter and characterised 
by its properties. Numerical simulations routinely find that the 
dark-matter distribution in galaxy clusters is expected to exhibit 
universal properties, for example its radial density profile and its 
degree of substructure. For our understanding of the nature of 
dark matter, it is important to test whether the mass distribution 
in real clusters confirms the expectations from simulations. 
Galaxy clusters are much less affected by baryonic physics than 
galaxies since the baryonic cooling time exceeds their lifetime 
except in their cores. Thus, they represent an important probe 
for the nature of dark matter and play a key role in testing our 
current understanding of cosmic structure formation. 
A growing number of increasingly wide surveys in a broad range 
of wavebands provide or will soon provide precise information 
on large galaxy-cluster samples. For instance, one of the goals of 
the ongoing Cluster Lensing And Supernova survey with Hubble 
(CLASH, (Postman et al. 2012)) is the mapping of dark matter 
in clusters based on strong and weak gravitational lensing. The 
Dark Energy Stirvey (DBS, (Frieman & Collaboration 2013)), 
having started in September 2012, will combine several probes 
of Dark Energy, and ESA's Euclid mission (Laureijs et al. 201 1) 
will mainly focus on weak gravitational lensing measurements 
and galaxy clustering, but will also provide data on galaxy clus- 
ters and the integrated Sachs-Wolfe effect. The Kilo Degree 
Survey (KiDS) aims to map the matter distribution in the uni- 
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verse by means of weak gravitational lensing and photometric 
redshift measurements (de Jong et al. 2013). 

Strong and weak gravitational lensing are widely used as 
an effective tool for the reconstruction of the projected 
mass or gravitational-potential distributions in galaxy clusters 
(Bartelmann et al. 1996; Cacciato et al. 2006; Merten et al. 2009; 
Coe et al. 2012; Brada et al. 2005; Bradac 2010; Meneghetti et al. 
2010; Medezinski et al. 2013). 

Lensing effects are due to light deflection only and thus (largely) 
insensitive to equilibrium and stability assumptions. They can be 
completely characterised by the scaled and projected Newtonian 
gravitational potential of the lensing matter distributions and 
thus most directly constrain the projected, two-dimensional 
gravitational potential. Non-parametric, adaptive methods have 
been developed and are now routinely being applied to re- 
cover cluster potentials. However, clusters provide a multitude of 
other observables through X-ray emission, the thermal Sunyaev- 
Zel'dovich effect and galaxy kinematics. In the series of articles 
which this paper belongs to, we are studying how all clusters 
observables could possibly be used to jointly constrain the pro- 
jected gravitational potential of galaxy clusters with as little prej- 
udice as possible. Combining all observables has the substan- 
tial advantage that all available information is bundled in single 
models, and that a range of linear scales covering approximately 
two orders of magnitude can faithfully be bridged. 

In this paper, we extend our previous studies towards including 
information from galaxy kinematics. The observables here are 
the line-of-sight velocity dispersions of the cluster galaxies as a 
function of cluster-centric radius. The relation between three- 
dimensional galaxy velocity dispersions and the gravitational 
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potential is governed in equilibrium by the Jeans equation. It re- 
sembles the equation of hydrostatic equilibrium, but contains an 
additional term taking account of a possible anisotropy in the ve- 
locity distribution. We describe the galaxy velocity dispersions 
here as an effective, possibly anisotropic pressure related to the 
matter density by an effective polytropic relation. Under this as- 
sumption, which we find well satisfied in simulations, the effec- 
tive galaxy pressure is related to the gravitational potential by 
a Volterra integral equation of the second kind, which can be 
solved iteratively. We can then proceed as in the previous papers 
of this series: We apply the Richardson-Lucy deprojection to the 
observable line-of-sight velocity dispersions to obtain a three- 
dimensional effective galaxy-pressure profile. This is then con- 
verted to the three-dimensional potential, from which the two- 
dimensional potential is found by projection. 
Our paper is organized as follows: We review in Sect. 2 the 
basic relations underlying our analysis, in particular our treat- 
ment of the Jeans equation and the Richardson-Lucy deprojec- 
tion. Numerical tests of our method based on kinematic data of 
a numerically simulated cluster and intended applications to ob- 
servational data are described in Sect. 3. We outline our conclu- 
sions in Sect. 4. 



2. Recovering the projected gravitational potential 
from the projected velocity dispersions 

2.1. Basic relations 

In order to incorporate measurements of the kinematics of clus- 
ter galaxies into reconstructions of the two-dimensional gravi- 
tational potential, we first require a relation between galaxy ve- 
locity dispersions and the three-dimensional gravitational poten- 
tial. The velocity dispersions are generally defined (see (Binney 
& Tremaine 1987; Schneider 2006)) as the mean squared devi- 
ations of the velocities of the cluster members from the mean 
velocity (v,) of the population: 



< = <v?) - {Vif 



(1) 



Measured velocity dispersions are density-weighted projections 
of the three-dimensional velocity dispersions along lines-of- 
sight through the cluster Throughout the paper, we orient the 
coordinate system such that the z-axis coincides with the line- 
of-sight. A projected velocity dispersion profile is constructed 
by averaging over concentric cylinderical shells drilled around 
the line-of-sight as a symmetry axis. This profile represents the 
input of our method. 

The final output of our algorithm will be the projected 
Newtonian potential of the lens, defined by 
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where (p is the three-dimensional gravitational potential of the 
cluster at distance D^ from the observer. The projected potential 
is given as a function of the two-dimensional position coordinate 
6 on the sky. 

The geometry of the problem is sketched in Fig. 2.1 (see also 
(Binney & Tremaine 1987)). For simplicity during the construc- 
tion of our method, we adopt a spherically-symmetric cluster 
model. All equations derived in the following will thus be for- 
mulated in spherical coordinates. 

The possible anisotropy of the velocity distribution will be de- 
scribed by the conventional anisotropy parameter /3 (Binney & 
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Fig. 2.1. Cluster geometry (see (Binney & Tremaine 1987)). 



Tremaine 1987), 



P^l- 
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Our algorithm consists of three main steps: 



- We deproject the observable, i.e. the velocity dispersions 
along the line-of-sight, into a deprojected quantity which, 
for reasons to be clarified later, will be called the effective 
galaxy pressure P. 

- Since this effective pressure P is related to the gravitational 
potential (p, we can solve the Jeans equation using symmetry 
assumptions and a formal analogy to a polytropic gas strati- 
fication. We thus obtain a relation between the gravitational 
potential and the effective galaxy pressure. This equation 
is a Volterra integral equation of the second kind for 0, which 
can be solved iteratively. 

- Having obtained the gravitational potential, we simply 
project it to find the two-dimensional potential ^. 

We will describe our deprojection method in Subsect. 2.2. The 
relevant three-dimensional quantities are the density p, the mean 
radial velocity dispersion cr]. and the gravitational potential (p. 
They are related via the Poisson and Jeans equations. In spherical 
symmetry, the Poisson and Jeans equations are 
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respectively. We note here that the second term on the left-hand 
side of Eq. (5) is the only formal difference to the hydrostatic 
equation for the hot intracluster gas. This difference is important 
because it complicates the solution of the Jeans equation consid- 
erably. 

For solving the Jeans equation, we formally identify the density- 
weighted radial velocity dispersion pcr^ with an effective galaxy 
pressure P and assume that it is related to the density by a poly- 
tropic relation. This assumption can be justified by the following 
calculation. 
A given radial density profile pir) implies the mass profile 



M{r) = 47r 
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and, by Poisson's equation, the gravitational-potential gradient 
50(r) GM{r) 



dr 



(7) 
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With this expression, the Jeans equation has the solution 
P(r) = Poexp(- r ^dxj 



X) 



dx\, 



(8) 



where the boundary conditions are set by the pressure Pq at the 
fiducial radius ro. The fiducial pressure Pq can be related to a 
fiducial density po via the virial theorem. 
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This way, the effective pressure profile of the cluster galaxies 
can be expressed in terms of the density and anisotropy profiles. 
We have evaluated Eq. (8) for two different mass density pro- 
files, viz. the Hernquist (Hernquist 1990) and NFW (Navarro 
et al. 1996) profiles, adopting the anisotropy profile derived by 
(Hansen & Moore 2006) from detailed numerical studies. As 
Fig. 2.2 shows, the effective pressure-density relation is nearly 
polytropic, i.e. an approximate power law. Thus, our assumption 
of an effectively polytropic galaxy stratification. 



P = P(; 



(10) 



seems appropriate. 







Fig. 2.2. The relation (8) between the effective pressure and the density 
is show here for the Hernquist and NFW density profiles, adopting the 
anisotropy parameter proposed by (Hansen & Moore 2006). The very 
nearly straight lines (note the logarithmic axes) demonstrate that the 
assumption of a polytropic relation is well justified. 



With this result, we return to the Jeans equation (5), where 
we express the density by the effective pressure by means of 
Eq. (10). For brevity, we further abbreviate 
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to cast the Jeans equation (Eq. (5)) into the form 
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This linear, inhomogeneous, first-order differential equation 
with non-constant coefficients can straightforwardly be solved, 
e.g. by variation of constants. The solution 



p^ - -(f{r) + exp 



-2r^dx 

•Jr„ X 



+ 2 



I dy— (^(y)expl2 I — dxj 



(13) 



is the unique relation between the effective pressure and the 
gravitational potential we were aiming at. 
Equation (13) is a Volterra integral equation of second kind 
(Abramowitz & Stegun 1972) that can be solved iteratively. In 
this way, we find a relation between the effective pressure p and 
the gravitational potential ip. Projection along the line-of-sight 
leads to the two-dimensional potential i/'. 

2.2. Deprojection 

Observations yield the line-of-sight projection of the galaxy ve- 
locity dispersion only. Three-dimensional velocity dispersions 
can be obtained by deprojection, which can be achieved e.g. by 
Richardson-Lucy deconvolution (Richardson 1972; Lucy 1974). 
As detailed in (Konrad et al. 2013), a line-of-sight-projection 
g(s) is assumed to be related to a three-dimensional quantity /(r) 
by 

g(s) = Jdz/( V^2 + ,2j ^ jdrf(r)K(s\r) (14) 

with a projection kernel K{s\r). In spherical symmetry, with the 
anisotropy parameter jS(r), the projection kernel for the velocity 
dispersion is 
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For isotropic velocity dispersions or an isotropic gas pressure, 
the final factor in parentheses in Eq. (15) is unity. The kernel in 
Eq. (15) can be easily derived from Fig. 2.1 as shown by (Binney 
& Tremaine 1987). 

Provided that g(s), fir) and K{s\r) are normalised', the follow- 
ing iterative scheme follows from Bayes' theorem (see (Lucy 
1974; Konrad etal. 2013)): 
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Thus, starting from a suitably guessed, normalised function 
fo(f), the scheme given by Eqs. (16) and (17) allows to recover 
the three-dimensional function f{r) from its two-dimensional 
projection g{s), assuming the symmetry incorporated into the 
projection kernel K{s\r). Experience shows that even for guess 
functions /o(r) wildly different from the true f(r), the method 
converges surprisingly quickly. 

For applying this approach to reahstic observational data con- 
taining small scale fluctuations due to background or instrumen- 
tal noise, a regularisation term should be introduced as discussed 



The kernel K(s\r) is normalised integrating over s. 
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in (Lucy 1994). It was shown there that the deprqjection algo- 
rithm can be interpreted as the resuh of a maximum-HkeUhood 
problem obtained by variation, 



fM{r)^fi{r) + Mr) 



^Him r. ...^nm 



fi(r) J 



drfiir)- 



Sfiir) J ■' Sfiir) 
where the functional derivatives of the functional H[f] 
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occur. 

Regularisation can be introduced by adding a functional S [/], 



H[f]^Q[f]^H[f]+aS[f] 
which can be of entropic form 
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Here, xif") is a smooth prior function chosen to suppress small 
scale fluctuations. A suitable prior can be the smoothed version 
of the deprojection result from the previous iteration step. This 
choice, known as the floating default (see (Home 1985; Lucy 
1994)), is defined by 



X(r) 



^ J dr' P(r\r')f(r') , 



(22) 



with a normalised, usually sharply peaked convolution kernel 
P(r\r') symmetric in r - r'. We use a properly normalised 
Gaussian with a smoothing scale L, 
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Replacing the variation of H[f] in Eq. (18) by the variation of 
Q[f] from Eq. (20), we obtain 



fi^i(r) = fi(r) 
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which we shall use henceforth. Since we are working with dis- 
cretised data sets, the integrals in Eq. (24) need to be approxi- 
mated by sums. 



3. Numerical tests 

We now proceed to demonstrate that it is possible to recover 
the projected gravitational potential i/' of a galaxy cluster from 
the measured velocity dispersions projected along the line-of- 
sight. The algorithm described in Subsect. 2.1 assumes spherical 
symmetry and a polytropic relation between the effective galaxy 
pressure and the density. 

Feeding the projected velocity dispersions into the Richardson- 
Lucy deprojection, we obtain the effective pressure P related to 
the gravitational potential (p by Eq. (13). Once this is achieved, 
the gravitational potential just needs to be projected along the 
line-of-sight. 



3.1. The data 

For testing our algorithm with simulated data, we require a 
velocity-dispersion profile projected along the line-of-sight and 
a two-dimensional gravitational potential obtained indepen- 
dently. We obtain such data from one of the clusters (denoted gl) 
in the hydrodynamically simulated sample described in (Saro 
et al. 2006) and used previously in (Meneghetti et al. 2010). 
We start from a catalogue listing the Cartesian coordinates and 
the three Cartesian velocity components of simulation particles 
tracking the motion of cluster galaxies. All information neces- 
sary for the kinematic analysis described in Subsect. 3.2 can be 
extracted from this catalogue. 

In addition, we converted the surface-mass density of the clus- 
ter into the two-dimensional potential, solving Poisson's equa- 
tion by means of a fast-Fourier transform. The surface-density 
map and the two-dimensional gravitational potential are shown 
in Fig. 3.1. Evidently, the potential is considerably smoother 
than the mass density. 

To enforce axial symmetry, we azimuthally average the gravita- 
tional potential shown in Fig. 3.1 (b) around the centre, chosen 
to be the point with the deepest (most negative) potential. In or- 
der to make it comparable to the outcome of our algorithm, it 
can be passed into the RL deprojection (with/? = 0) yielding the 
gravitational potential which can then be appropriately shifted 
and normalised. 



3.2. Testing the algorithm 

Assuming spherical symmetry, we transform from Cartesian to 
spherical polar coordinates using 

r = Jx^ + y^ + z^ , 
X y z 

Vr = V^- + Vy- + Vj- , 

r r r 

y X 
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" y x^ H- y^ r ^Jx^ + y^ 
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The number-density profile of cluster galaxies is obtained within 
radial bins r, centred on the cluster centre chosen to contain 
equal numbers of galaxies. The number counts are then con- 
verted to a number-density profile by weighting with the inverse 
volume of each radial shell. 



n(r) oc 



counts 



(26) 



Figure 3.3 reveals that the number-density profile for this partic- 
ular cluster essentially follows a power law. In a first approxima- 
tion sufficient for our purposes, galaxy biasing and variations of 
the galaxy mass function are neglected, allowing us to adopt the 
number density as a direct tracer of the mass density. 
Given the number-density profile and the radial bins, a mean 
galaxy velocity can be obtained in each radial shell. The vari- 
ance of the velocities about this mean yields the galaxy velocity- 
dispersion profiles cr^, (A and cr^ in the radial, azimuthal and 
polar directions. These quantities enable us to derive the profile 
of the anisotropy parameter according to the definition in Eq. (3). 
The result is shown in Fig. 3.2. Despite the substantial scatter, we 
fit the two-parameter model 



/3(r)^/3c. 
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(a) 




Fig. 3.1. (a) Surface-mass density of the simulated cluster gl in the x-y plane, (b) Two-dimensional gravitational potential after obtained from the 
surface-density map solving Poisson's equation via fast-Fourier transform. 



proposed by (Mamon & Lokas 2005), generalized by (Tiret et al. 
2007) and used in (Mamon et al. 2013). This fitted profile will 
be used as the mean anisotropy-parameter profile in the recon- 
struction, despite its being a rough approximation. 
We weigh the radial velocity dispersions cr^ with the number- 
density and obtain the effective galaxy pressure. Projecting this 
quantity along the line-of-sight. 



-r 



dr p,{r)K(s\r) , 



(28) 



with the projection kernel K(s\r) defined in Eq. (15), we finally 
complete the input for our method which would be a proper ob- 
servable provided by observations. Note that the observed ve- 
locity dispersions are quantities projected along the line-of-sight 
and are thus implicitly weighted by the number density of galax- 
ies along the line-of-sight. 

Figure 3.4 is a double-logarithmic plot of the effective galaxy 
pressure vs. the density, confirming that our polytropic assump- 
tion is reasonable. The mean polytropic index y, introduced in 
Sect. 2.1, is estimated from it by linear regression. 
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Fig. 3.3. The number density profile of the simulated cluster galaxies is 
show vs. the radius. It is obtained by galaxy number counts in spherical 
shells. 
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Fig. 3.2. Radial profile of the anisotropy parameter fi(r) defined in 
Eq. (3), obtained from the simulated cluster data. The two-parameter 
model by (Mamon & Lokas 2005) is fitted to the data points. We find 
jSoo = 0.40±0.5 and r^ = 0.08±0.74. Concerning this result, a vanishing 
P would also be allowed. 



Figure 3.5 illustrates our complete algorithm. In the top-left 
panel, we show the normalised, line-of-sight projected ve- 
locity dispersions as a function of the projected radius and 
weighted with the galaxy number density. In the top-right panel, 
the normalised effective pressure profile obtained from the 
Richardson-Lucy deprojection algorithm is displayed. We then 
invert Eq. (13) to obtain the three-dimensional, Newtonian grav- 
itational potential shown in the bottom-left panel. In the bottom- 
right panel, we finally show the derived two-dimensional gravi- 
tational potential. 

Figure 3.5 provides an overview over the first test we performed 
on our algorithm and shows all steps needed to reconstruct the 
two-dimensional potential from the line-of-sight projected ve- 
locity dispersions. A comparison between our reconstructed pro- 
files and the true three- and two-dimensional gravitational poten- 
tials is now in order. We define as "true" the potentials extracted 
directly from the numerically-simulated cluster 
Figure 3.6 (a) shows the comparison between the reconstructed 
and the true gravitational potentials, while Fig. 3.6 (b) compares 
the reconstructed and the true projected potentials. In both pan- 
els, the blue curve shows the true potential extracted directly 
from the simulation, while the red curve shows our reconstruc- 
tion results. For easier comparison of both potentials, their zero 
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Fig. 3.5. Results of the four different steps comprising our algorithm, (a) The input to our pipeline is the normalised, line-of-sight projected velocity- 
dispersion profile as a function of the projected radius and weighted by the galaxy number density, (b) The Richardson-Lucy deprojection algorithm 
yields an effective galaxy-pressure profile, (c) Solving the Volterra integral equation (Eq. (13)), we obtain the three-dimensional, Newtonian 
potential, (d) The last step consists in projecting the gravitational potential of panel (c) to find the two-dimensional potential. 



points and normalisations are adjusted^. To be consistent with 
(Konrad et al. 2013) and (Majer et al. 2013), we decided to nor- 
malise both functions to unity and fix (pir^ax) - 0. In Fig. 3.7, 
the relative deviation 



Ai//(s) 



|lA™(i)| 



(29) 



between the reconstructed and the true two-dimensional poten- 
tials is shown as a function of projected radius. The deviation 
strongly increases at large radii, although it remains below 10 % 
within 1.2/i^'Mpc, which is similar to the virial radius of the 
cluster. 



^ Note that due to the normalisation constraint of the Richardson- 
Lucy deprojection, we obtain a function proportional to 0. The poten- 
tial extracted from the simulation has to be normalised in the same way. 
Since both potentials are obtained from different datasets providing dif- 
ferent spatial boundaries, they need to be appropriately shifted (gauged) 
by adding a constant factor. 



3.3. Application to observed data 

The numerical tests presented above were run on line-of-sight 
projected data extracted from a simulated galaxy cluster Starting 
from spatial Cartesian coordinates and velocity components of a 
sample of member galaxies, the effective galaxy pressure was 
obtained by weighting the squared radial velocity dispersion 
with the number density. Afterwards, the required observable 
was constructed by projection along the line-of-sight. The ef- 
fective polytropic index -y and the profile /3(r) of the anisotropy 
parameter could be easily obtained from the simulated sample. 
Unfortunately, none of the latter two quantities is directly ac- 
cessible to observations. Therefore a method needs to be pro- 
posed that allows to model the effective polytropic index and the 
anisotropy profile. 

Major efforts have been undertaken in recent years to either con- 
strain the anisotropy profile from observational and simulated 
data (Lemze et al. 2009; Mamon et al. 2013) or to identify a 
general relation describing it (Mamon & Lokas 2005; Hansen 
& Moore 2006). The main objective of our work is the recon- 
struction of the projected gravitational potential of a galaxy clus- 
ter. However, we can envisage an inverse application of our al- 
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Fig. 3.6. The reconstructed gravitational potentials in (a) three and (b) two dimensions are plotted as functions of radius and compared with the 
true potentials. In each panel, the blue points show the true potential determined from the convergence map, while the red points show the result 
of our reconstruction method. The corresponding relative deviation as a function of radius is shown in Fig. 3.7. 
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gorithm. By an operation analogous to (Lemze et al. 2009), it 
is possible to constrain /3(r) using simultaneously information 
from gravitational lensing and galaxy kinematics. The existence 
of multiple constraints on the two-dimensional potential iff also 
allows developing an iterative method for constraining /3{r) and 

r- 

Within the CLASH survey project, considerable progress has 
been made in reproducing the observed CLASH clusters by sim- 
ulations (Meneghetti et al. 2013, in prep.). These results would 
possibly enable us to guess a polytropic index and an anisotropy 
profile which could then be combined with kinematic data from 
observations in order to improve our reconstruction. 



4. Conclusions 

With the study presented here, we have continued our series 
of papers aiming at parameter-free reconstructions of projected 
gravitational potential of galaxy clusters. The essential goal of 



our approach is thus to construct methods to relate all cluster ob- 
servables to the two-dimensional gravitational potential. In two 
preceding papers, we have shown how this can be achieved with 
the two cluster observables based on the hot intracluster gas, i.e. 
X-ray emission and the thermal Sunyaev-Zel'dovich effect. In 
this paper, we address the problem how the observed line-of- 
sight projected velocity dispersions of cluster galaxies can be 
converted to the projected gravitational potential. We propose a 
solution which follows closely our interpretation of the observ- 
ables provided by the intracluster plasma, but necessarily differs 
from it in the crucial aspect that the velocity dispersion of the 
cluster galaxies, unlike the gas pressure, can be anisotropic. 
The algorithm we propose rests on the following crucial assump- 
tions. First, we assume that the effective galaxy pressure, by 
which we mean the product of the galaxy number density and 
the radial velocity dispersion, can be related to the density it- 
self by a polytropic relation, i.e. a power law. We test this as- 
sumption with different density profiles and functional forms 
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of the anisotropy parameter and find it well satisfied. Second, 
adopting this polytropic relation, we solve the radial component 
of the Jeans equation, relating the effective galaxy pressure to 
the potential gradient. This solution can be analytically given in 
the form of a Volterra integral equation of the first kind for the 
three-dimensional gravitational potential, which can be solved 
by iteration. Thus, we establish a relation between the three- 
dimensional gravitational potential and the eff'ective galaxy pres- 
sure. 

The effective galaxy pressure itself can be obtained from the 
observable, line-of-sight projected galaxy velocity dispersions 
by Richardson-Lucy deprojection. The two-dimensional gravi- 
tational potential can finally be found by straightforward projec- 
tion along the line-of-sight. 

We have tested this algorithm on a simulated galaxy cluster in 
which galaxies have been identified. This cluster is part of a 
sample of numerical hydrodynamical simulations described in 
(Saro et al. 2006) and used in (Meneghetti et al. 2010). The 
anisotropy profile [i{r), as well as the mean effective polytropic 
index y, were obtained directly from the simulation. Although 
these quantities are hardly directly measurable from observa- 
tions and it was up to now impossible to give general prescrip- 
tions of their behaviour, it may well be justified to calibrate them 
on numerical simulations without introducing an unacceptable 
bias. 

The numerical simulation shows that the anisotropy profile 
shown in Fig. 3.2 may be poorly constrained by the given simu- 
lation. The fit performed by the model suggested by (Mamon 
& Lokas 2005) covers a wide range of possible behaviour. 
Nonetheless, reconstructing the two-dimensional potential based 
on the nominal fit yields a result that closely follows the numeri- 
cal expectation out to « 1 .2 /j^'Mpc, as demonstrated by the rel- 
ative deviation of the reconstructed from the true projected grav- 
itational potential (see Fig. 3.7). The deviation remains below 
10 % in this radial range. Figure 3.4 confirms the approximate 
treatment of our system as a polytropic stratification. 
The main limitation of our algorithm so far lies in the assumption 
of spherical symmetry, which we plan to relax by extending our 
method to the axially (rather than spherically) symmetric case. 
We are planning to proceed by applying the algorithm devised 
here to observational data that could possibly be provided by the 
CLASH survey. In addition, we will focus on unifying the three 
algorithms presented in (Konrad et al. 2013; Majer et al. 2013) 
and this paper such as to allow a single, joint reconstruction of 
the two-dimensional potential of galaxy clusters using lensing. 
X-ray, Sunyaev-Zel'dovich and kinematic data. 
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